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Minute  amounts  of  ampiphillic  molecules  at  an  air-water  inter¬ 
face  can  have  striking  effects  on  both  the  damping  and  generation  of 
capillary  waves.  The  theory  of  these  effects  within  linearized 
hydrodynamics  is  reviewed.  The  conventional  wisdom  is  that  oil  calms 
troubled  waters  by  increasing  the  rate  of  capillary  wave  damping. 
Although  surface-active  films  can  increase  the  damping  rate  by 
factors  of  two  to  four,  a  potentially  more  important  effect  is  the 
suppression  of  capillary  wave  generation  by  the  wind.  In  particular, 
generation  by  the  components  of  fluctuating  Reynold's  stresses  in  the 
atmosphere  parallel  to  the  surface  can  be  reduced  by  an  order  of 
magnitude  or  more.  The  theory  predicts  a  more  modest  suppression  of 
the  effect  of  Reynold's  stresses  normal  to  the  surface. 
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1.0  INTRODUCTION 


The  calming  effect  of  oily  materials  on  water  waves  has  fasci¬ 
nated  scientists  and  seafarers  since  antiquity.  In  a  famous  experi¬ 
ment,  Benjamin  Franklin  observed  that  a  teaspoon  of  oil  spread 
rapidly  to  cover  approximately  half  an  acre  on  a  pond,  producing  an 
area  "smooth  as  a  looking  glass. Assuming  that  the  teaspoon  con¬ 
tained  about  5  cm^  of  oil,  we  can  calculate  that  Franklin's  film  was 
about  20  Angstroms  (2  x  1 0-^  cm)  thick.  Recently,  sun-glint  photo¬ 
graphs  of  the  western  Mediterranean  Ocean  taken  from  the  Space 
Shuttle  have  revealed  a  complicated  pattern  of  quiescent  regions, 
which  reflect  the  sun  specularly,  as  well  as  regions  of  diffuse 

p 

reflection,  presumably  due  to  a  profusion  of  capillary  waves.  The 
regions  of  specular  reflection  have  been  attributed  to  sea  slicks, 
possibly  of  natural  origin.  (Diatoms,  for  example,  contain  droplets 
of  oil  to  assist  in  flotation  and/or  as  an  emergency  food  supply.-^) 

As  representative  surface-active  materials,  we  may  consider 
saturated  fatty  acids  (lipids)  with  the  chemical  formula 
CH^  -  (CH2)n  ~  COOHi  where  n  is  usually  between  10  and  20.  The  polar 
head  group  is,  for  all  practical  purposes,  infinitely  soluble  in 
water;  the  hydrocarbon  tail,  however,  is  hydrophobic,  because  it 
disrupts  the  delicate  network  of  hydrogen  bonds  which  characterize 
short  range  order  in  liquid  H20.  The  competing  hydrophillic  and 
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hydrophobic  tendencies  tend  to  pin  such  "ampiphill ic"  molecules 

normal  to  an  air-water  interface,  with  the  tails  directed  towards  the 

0 

atmosphere.  Since  ampiphiles  typically  measure  20  A  along  their  long 
axis,  Franklin's  film  was  one  molecular  layer  thick. 

How  can  such  a  thin  layer  cause  a  pronounced  change  in  the  wave 
properties?  The  well-known  reduction  in  surface  tension  due  to  the 
surfactant  produces  only  a  minor  change  in  the  velocity  of  the 
waves.  Hydrodynamic  theories  of  how  surface-active  films  affect 
capi 1 lary-gravity  wave  damping  have  been  worked  out  by  Reynolds,^ 
Levich,^  and  Dorrestein.^  who  find  effects  which  are  most  pronounced 
in  the  capillary  wave  part  of  the  spectrum.  The  basic  idea  is  that 
the  finite  compressibility  of  a  lipid  monolayer  modifies  the  boundary 
conditions  at  an  air-wate  ’  interface.  This  change  in  boundary  condi¬ 
tions  will  affect  the  flow  below  the  surface  down  to  a  viscous 

7 

penetration  depth  5, 

6  =  /2v/w  ,  (1.1) 

where  v  is  the  kinematic  viscosity  of  water  and  w  is  the  angular 
frequency  of  the  wave.  Since  5  is  of  order  10  -  1 0  -1  cm  for  capil¬ 

lary  waves  with  wavelengths  of  order  .1  to  1  cm,  surfactant  molecules 
influence  the  flow  down  to  a  depth  which  vastly  exceeds  their  3ize. 
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When  surfactants  are  present,  the  vorticity  in  the  viscous  penetra¬ 
tion  depth  increases,  increasing  the  damping  by  factors  as  large  as  2 
to  4  relative  to  a  clean  surface. 

Although  remarkable,  it  seems  unlikely  that  this  enhanced  damp¬ 
ing  could,  by  itself,  account  for  interfaces  which  are  "as  smooth  as 
a  looking  glass"  when  covered  with  an  oily  film.  On  a  clean  surface, 
a  0.5  cm  capillary  wave,  for  example,  decays  in  a  time  which  is  more 
than  twenty  times  its  basic  period.  Capillary  waves,  once  generated, 
will  persist  for  several  wave  periods  even  when  the  damping  is 
enhanced  by  surface-active  films.  The  damping  can  never  become  very 
large,  because  the  dissipation  is  confined  to  the  viscous  penetration 
depth,  which  depends  only  weakly  on  the  surfactant  concentration. 

The  kinetic  energy  of  the  flow  is  distributed  over  a  much  larger 
volume,  whose  depth  is  comparable  to  the  wavelength.  The  decay  rate 
is  increased  by  a  most  a  factor  three  even  in  the  extreme  case  of 
films  which  are  incompressible  or  have  infinite  viscosities. 

Sea  slicks  can  also  inhibit  the  generation  of  capillary  waves  by 

O 

the  wind.  For  clean  surfaces.  Lamb  discusses  a  Kelvin-Helmholtz- 
like  instability  leading  to  wave  growth  above  a  critical  wind  speed 
of  6.5  m/sec. ^  Lamb  points  out,  however,  that  waves  are  apparently 
generated  by  other  means  at  considerably  lower  velocities.  A  number 
of  alternative  theoretical  approaches  are  available.10 


Using  one 


such  approach,  Miles  finds  a  minimum  wind  speed  of  about  1  m/sec  for 
a  clean  surface,  and  estimates  that  an  incompressible  film  increases 
this  threshold  by  an  order  of  magnitude.11 

In  this  report,  we  study  a  different  mechanism  for  wave  genera- 
tion,  namely  turbulent  fluctuations  of  the  wind  about  its  mean 
velocity.  If  a  wave  with  a  particular  wavevector  and  frequency  is 
already  present,  it  will  be  acted  upon  most  strongly  by  those  compo¬ 
nents  of  the  fluctuating  atmospheric  Reynold's  stress  with  matching 
wavevectors  and  frequencies.  This  particular  generation  mechanism 

Q  1  o 

has  been  explored  for  clean  surfaces  by  Lamb7  and  by  Phillips. 

3y  calculating  the  response  of  capillary  waves  to  fluctuating 
Reynold's  stresses  within  linearized  hydrodynamics,  we  show  that 
surfactants  produce  a  pronounced  reduction  in  the  response  to 
stresses  parallel  to  the  surface.  If  only  tangential  stresses  are 
present,  a  modest  amount  of  surfactant  can  reduce  the  root  mean 
square  amplitude  of  capillary  waves  by  a  factor  of  10  or  more  rela¬ 
tive  to  a  clean  surface;  the  root  mean  square  amplitude  tends  to  zero 
ir.  toe  limit  of  an  incompressible  film.  There  is  also  a  reduction  in 
tr.e  response  to  fluctuating  normal  stresses,  although  here  the  effect 
is  only  a  f actor  of  2  or  3-  The  different  responses  to  parallel  and 
normal  stresses  agrees  with  a  suggestion  by  Dorrestein.^  Little  is 
-mown  about  the  -dative  strengths  of  fluctuating  tangential  and 


normal  stresses  near  the  ocean  surface.  If  these  stresses  are  equal, 
the  root  mean  square  wave  amplitudes  will  be  reduced  by  a  factor  of 
about  6.  If  tangential  stresses  dominate,  however,  the  suppression 
of  capillary  waves  by  surfactants  could  be  much  more  effective, 
especially  for  films  with  low  compressibilities. 

In  Section  2.0,  we  review  some  elementary  properties  of  lipid 
films  which  will  be  needed  later  in  the  report.  The  hydrodynamic 
theory  of  damping  of  capillary-gravity  waves  by  surfactants  is 
reviewed  in  Section  3.0.  This  theory  allows  the  damping  to  be  pre¬ 
dicted  directly  from  laboratory  measurements  of  pressure-area 
isotherms;  in  some  cases,  knowledge  of  film  viscosities  may  also  be 
necessary.  Finally,  in  Section  4.0,  we  calculate  the  effect  of  lipid 
monolayers  on  the  response  to  fluctuating  atmospheric  Reynold's 
stresses  at  the  sea  surface. 


2.0  PROPERITIES  OF  SURFACE-ACTIVE  FILMS 


Figure  1  shows  a  typical  experimental  setup  for  studying 
surface-active  films  in  the  laboratory.13  Lipid  moelcules  are  con¬ 
fined  to  a  region  of  area  A  by  a  float.  The  air-water  surface 
tension  oH  Q  «  71*  dyne/cm  is  reduced  by  the  presence  of  the  lipid 
monolayer.  As  a  result,  the  film  exerts  a  pressure  n  on  the  float  as 
it  tries  to  spread  and  cover  the  clean  water  surface.  Upon  measuring 
the  pressure  necessary  to  confine  the  film  to  different  areas  at 
constant  temperature  T,  one  obtains  a  pressure-area  isotherm  like 
that  shown  in  Figure  2.  It  is  easily  shown  that  the  pressure  tt  gives 
the  reduction  in  the  surface  tension  due  to  the  presence  of  the 
film,13 


o  = 


! 


(2.1  ) 


where  o  is  the  surface  tension  with  the  film  present.  To  a  first 
approximation,  pressure-area  isotherms  at  temperature  T  are  often 
described  by  a  two-dimensional  ideal  gas  equation  of  state, 


where  N  is  the  total  number  of  lipid  molecules  in  the  film.  In  real 
films,  the  pressure  rises  more  steeply  with  decreasing  area  than 
predicted  by  Equation  (2.2),  due  to  hard  core  repulsions  between 
lipids. 

At  high  lipid  densities,  there  can  be  anomalies  in  pressure-area 
isotherms  due  to  phase  transitions  in  the  film.  Many  different 
phases  are  possible,  including  solids,  liquids,  gases,  tilted 
liquids,  expanded  liquids,  tilted  solids,  etc.1^  In  principle,  there 
could  be  enhanced  damping  of  surface  waves  if  the  periodic  pressure 
variations  in  a  wave  were  to  drive  a  film  back  and  forth  across  one 
of  these  transitions.  It  seems  improbable,  however,  that  naturally 
occuring  slicks  would  appear  exactly  at  a  density  corresponding  to 
one  of  these  phase  transitions.  The  densities  in  sea  slicks,  more¬ 
over,  are  probably  those  of  a  dilute  liquid  or  moderately  dense  gas,2 
for  which  most  of  these  transitions  are  absent. 

Besides  reducing  the  surface  tension,  surface-active  monolayers 
also  cause  the  air-water  interface  to  behave  like  an  elastic 
membrane.  When  the  surface  distorts,  there  are  now  elastic  restoring 
forces  in  the  plane  of  the  film.  The  elasticity  arises  because  the 
surface  tension  depends  on  the  concentration  r  =  N/A  of  surfactant. 
Consider,  for  example,  the  behavior  of  the  teaspoon  of  oil  discussed 
in  the  Introduction.  The  film  initially  spreads  rapidly  with  a 


velocity  of  about  10  cm/sec. ^  It  eventually  stops  spreading,  how¬ 
ever,  at  monolayer  thicknesses,  because  further  dilution  of  the  film 
would  increase  the  surface  tension. 

To  calculate  the  force  which  stops  the  spreading,  consider  the 
free  energy  of  a  film  of  area  A, 
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oA  + 


(2.3) 


where  A  is  the  area  which  is  free  of  film.  Let  us  compute  the  change 
in  free  energy  associated  with  a  small  change  in  the  film  area, 

A  -»  A  +  6A ,  A  -»  A  -  6A .  The  surfactant  concentration  changes 
according  to 
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Adding  up  the  changes  in  Equation  (2.3),  we  see  that  the  change  in 
free  energy  is 

6F  -  [o  -  oH  Q  -  r  6A  .  (2.6) 

Stability  of  the  film  requires  that  this  change  vanish,  which  means 
that  the  "spreading  pressure"  it  «  o„  n  -  o  must  be  balanced  by  the 

S 

term  arising  from  density  dependence  of  the  surface  tension, 
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(2.7) 


The  quantity 


B  3  '  r  2  °  '  (2.8) 

acts  like  an  elastic  bulk  modulus  of  the  film,  and  will  play  an 
important  role  in  the  hydrodynamic  theory  discussed  in  the  next  two 
sections.  Note  that  for  a  thick  film,  o  is  independent  of  f,  so  that 
this  bulk  modulus  vanishes.  For  a  film  which  has  spread  to  monolayer 
thicknesses  with  no  external  constraints,  Equation  (2.7)  shows  that  B 
equals  the  spreading  pressure.  Films  can  also  be  compressed,  by 

artificial  barriers  like  that  shown  in  Figure  1,  or  by  internal  waves 
15 
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at  sea. 


In  this  case,  we  have,  using  Equation  (2.1) 


Equation  (2.9)  shows  that  B  is  the  inverse  two-dimensional  compressi¬ 
bility  and  can  be  determined  by  differentiating  pressure-area 
isotherms.  For  an  ideal  gas  (which  would  spread  to  infinite  dilution 
unless  it  were  confined),  we  have 

B  =  kgT/a  ,  (2.10) 

where  a  =  A/N  is  the  area  per  molecule. 


3.0  DAMPING  OF  CAPILLARY-GRAVITY  WAVES  BY  SURFACTANTS 


Assume  the  mean  air-water  interface  is  normal  to  the  z-axis,  and 
consider  the  propagation  of  a  deep  water  surface  wave  in  the  x-direc- 
tion  (See  Figure  3).  Following  the  methods  of  References  ^-6,  we 
write  x-  and  z-components  of  the  velocity  as 


v  =  3  <(>  -  3  \|) 
X  X  2 


(3.1a) 


v  =  3  d>  +  3  ill 
y  z  z 


(3.1b) 


where  <j>  and  ip  are  functions  of  x,  z,  and  time  t.  The 
function  4>  describes  inviscid,  incompressible  potential  flow  and 
satisfies  Laplace's  equation. 


V^d,  =  0  . 


(3.2) 


The  linearized  incompressible  Navier-Stokes  equations, 


3  v  -  -  —  $ 

r  n 


2  + 
v  V  v 


will  be  satisfied,  provided  i|i  satisfies 


3t* 


v  V^i ip 


(3.H) 


where  v  is  the  kinematic  shear  viscosity  and  pQ  is  the  density. 
Restricting  our  attention  to  deep  water  waves,  we  write  41  and  ij;  as 
periodic  functions  of  x  and  t  which  are  exponentially  damped  in  the 
negative  z  direction. 


4>(x,  z,  t) 


.  ikx  +  kz 
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(3.5a) 


i|i(x,z,t) 


_  ikx  +  Zz  -  iuit 
Ce  , 


(3.5b) 


where  A  and  C  are  complex  amplitudes  which  will  be  determined  by 
boundary  conditions.  The  undisturbed  surface  is  taken  to  be  z  =  0. 
The  ansatz  (3.5a)  for  <p  automatically  satisfies  Laplace's  equation; 
ip  will  satisfy  Equation  (3-^)  provided  that  Z  is  related  to  k  by 


2  2 

l  «  k  -  iui/v  . 


(3-6) 


The  form  of  4  and  i|>  means  that  the  velocities  v  and  v  may  be 

x  y 


written 
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(3.7b) 


We  require  that  <j>  satisfy  the  inviscid  Navier-Stokes  equations, 
which  means  that  the  pressure  is  given  by. 


P  -  'P08z  +  iwpQAe 


ikx  +  kz  -  iu)t 


(3.8) 


where  g  is  the  gravitational  constant,  and  we  take  the  pressure  at 
the  undeformed  interface  to  be  zero.  Boundary  conditions  will  ulti¬ 
mately  be  specified  at  the  displaced  interface  z  =  c(x,t).  Since 
3.  C  =  v_  in  the  limit  of  small  velocities  and  displacements,  we  have 
using  Equation  (3.7b)  that 


C(x,t) 


kA+ikB  ikx  -  iuit 

- : -  e 

- 1(1) 


(3.9) 


Because  $  obeys  Laplace's  equation,  its  contribution  to  the 
viscous  damping  in  the  Navier-Stokes  equations  vanishes.  Dissipation 
comes  from  the  rotational  part  of  the  flow,  which  is  controlled 


by  Hi: 


We  shall  be  primarily  interested  in  waves  with  wavelengths  X  =  2n/k 


which  are  much  larger  than,  say,  10  ^  cm.  For  such  waves,  which  are 
only  weakly  damped,  it  is  an  excellent  approximation  to  neglect  the 
k2  term  in  Equation  (3.6)  and  write 


The  decay  of  the  function  i|i  in  the  z-direction  is  controlled  by  the 
real  part  of  l,  which  is  just  the  reciprocal  of  the  viscous  penetra¬ 
tion  depth  6  defined  in  the  introduction.  Equation  (3.11)  is  a  good 
approximation  for  Si  whenever  the  wavelength  is  much  greater  than  the 
viscous  penetration  depth, 


X  >>  6  . 


(3.12) 


Dissipation  occurs  only  within  the  viscous  penetration  depth;  its 
strength  is  given  by  the  complex  amplitude  C  in  Equation  (3.5),  which 
is  in  turn  determined  by  the  boundary  conditions.  The  flow  is 
guaranteed  to  be  potential  flow  at  depths  much  greater 
than  6,  regardless  of  the  boundary  conditions  (see  Figure  3). 


At  the  interface  there  are  boundary  conditions  on  both  the 
normal  and  tangential  forces.  The  normal  boundary  condition  is 
obtained  by  equating  the  zz-component  of  the  fluid  stress  tensor, 
evaluated  at  z  =  c(x,t),  to  the  usual  capillary  restoring  force, ^ 


(-P 


2n  3  v 
z  z 


>1. 


(3.13) 


In  the  linearized  theory,  surfactants  only  affect  the  normal  boundary 
condition  by  changing  the  surface  tension  o  according  to 
Equation  (2.1).  The  tangential  boundary  condition  follows  from 
requiring  that  the  xz-component  of  the  fluid  stress  tensor,  evaluated 
at  z  =  e(x,t),  match  the  forces  arising  from  the  flow  of  the  surfac¬ 
tant  within  the  surface.  In  analogy  with  three-dimensional 
compressible  fluids,  we  define  a  two-dimensional  viscous  stress 
tensor^ , 
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(3.14) 


where  i  and  j  denote  directions  in  the  (x,y)-plane,  and  and  i ; 
are  two-dimensional  shear  arid  bulk  viscosities.  We  assume  the  lipid 
molecules  are  constrained  to  move  at  velocities  given  by  the  bulk 
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expressions  (3.7)  evaluated  at  the  interface.  Since  vy  is  zero,  the 
force  due  to  the  surfactant  at  the  interface  is  directed  along  the 
x-axis  and  equals 
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(3.15) 


where  we  have  defined  a  composite  two-dimensional  viscosity  p^, 


(3.16) 


Upon  recognizing  that  it  is,  up  to  an  additive  constant,  just  the 
negative  of  the  surface  tension  (see  Equation  (2.1)),  we  see  that  the 
tangential  boundary  condition  becomes 


nO  v 


9  v  ) 
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(3.17) 


To  proceed  further,  we  need  to  know  the  space  and  time 
dependence  of  a  induced  by  the  wave  motion.  Assuming  that  o  changes 
by  virtue  of  changes  in  the  surfc  itant  concentration  F,  we  have 


so  we  need  to  know  r(x,t).  Assuming  that  r  undergoes  both  diffusion 
and  convection  at  the  interface,  we  haveJ 

3  r  +  V(vr)  =  d  v2r  .  (3.19) 

The  two-dimensional  diffusion  constant  D  wm  be  controlled  by  diffu¬ 
sion  of  the  lipid  polar  head  groups  at  the  interface  through  water. 

It  certainly  cannot  be  much  greater  than,  say,  ICT^  cm2/sec,  a  typi¬ 
cal  diffusion  constant  in  a  bulk  liquid.  Under  these  conditions,  the 
convective  part  of  Equation  (3.19)  dominates  over  diffusion.  We 
neglect  the  finite  solubility  of  lipid  molecules  in  water  which, 
according  to  Levich,-’  decreases  the  ability  of  surfactants  to  damp 
surface  waves.  Neglecting  D,  and  writing  r  in  terms  of  the 
deviation  r'(x,t)  from  its  equilibrium  value  r  , 

r(x,t)  =  rQ  +  r'(x.t)  ,  (3.20) 


we  have,  to  lowest  order  in  the  small  quantities  vx  and  f', 


or,  integrating  with  respect  to  time  and  differentiating  with  respect 
to  x, 


~ir°  a  2 
-  3  v 

u  xx 


(3.22) 


Combining  Equation  (3-22)  with  Equations  (3.18)  and  (3.17),  we 
see  that  the  tangential  boundary  condition  becomes 


n(3  v  +  3  v  )  I  =  (—  +  u_)  3  2  v 

z  z  xx'z=s  (o  2  x  x 


(3.23) 


where  we  have  introduced  the  surfactant  bulk  modulus  B  defined  by 
Equation  (2.8).  Note  that  the  bulk  modulus  can  be  regarded  as 
contributing  to  a  frequency- dependent  two-dimensional  viscosity.  The 


*  .  v  ,v  V“ 
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angular  frequency  of  a  .5  cm  capillary  wave  is  about  300  Hz.  For  a 

°2 

moderately  dense  film,  with  40  A  per  molecule,  the  ideal  gas  formula 
(2.10)  gives  B  =  10  dyne/cm.  Under  these  conditions,  thr  contribu¬ 
tion  of  the  bulk  modulus  to  the  effective  viscosity  is 


,'®C  . 


B  "2 

-'3x10  dyne-sec/cm 
w 


(3.24) 


Since  typical  intrinsic  film  viscosities  are  less  than  or  equal  to 
10*"4  dyne- sec/cm, 1 ^  we  shall  neglect  the  intrinsic  viscosity  in  most 
of  what  follows. 
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Setting  to  zero,  and  using  Equations  (3.7a),  (3.7b),  (3.8), 
and  (3.9)  in  the  boundary  conditions  (3.13)  and  (3.23),  we  find  two 
homogeneous  equations  for  the  complex  amplitudes  A  and  C,  namely 

(-w2  +  a)  2  -  2ivwk2)  A  +  (iui  2  +  2vk£u>)  C  -  0  ,  (3.25a) 

o  o 

2 

p  O  p  p  RV  0 

Uvuk*  +  iBk-3/po)  A  ♦  (u  ♦  2ivuk*  -  )  C  =  0  (3.25b) 


where 


“o<k) 


ok^/p 


(3.26) 


Simultaneous  solutions  are  only  possible  if  the  determinant  of 
coefficients  vanishes. 


A(k,w)  a  det 


2  2  2  2 
-(o  +  oj  -  2ivwk  ,  iu  +  2vk!,u) 

o  o 


pop  p  Rk  5. 

2vwk  +  2BkJ/p  ,  a)  +  2ivwk  -  - 

o’  P„ 


=  0 


(3.27) 


which  defines  implicitely  the  dispersion  relation  w  =  u(k)  of 
capillary-gravity  waves.  If  there  is  no  surfactant  present,  B  =  0, 


and  Equation  (3.27)  leads  in  the  limit  X  >>  6  to  the  standard 
result,^ 

u»(k)  *  (u0(k)  -  2ivk2.  (3.28) 

Suppose,  on  the  other  hand  B  =  °°,  so  that  the  film  is  incom¬ 
pressible.  This  is  the  limiting  case  discussed  by  Lamb.^  From 
Equation  (3.25b),  we  have  that 

ikA  -  5.B  =  0  .  (3.29) 

which,  in  view  of  Equation  (3.7a),  means  that  the  tangential 
component  of  velocity  vanishes,  vx  =  0.  As  shown,  for  example,  in 
the  book  by  Levich,^  the  solution  of  (3.27)  in  the  capillary  wave 
part  of  the  spectrum  (X  <  1.7  cm)  is  then 


where  we  have  neglected  a  small  shift  in  a  real  part  of  the 

frequency.  The  ratio  of  the  imaginary  part  of  (3.30)  to  the  damping 
2 

2vk  for  a  clean  surface  is 


Im  a) 


1 


(3.31  ) 


where  the  wavelength  in  the  last  line  is  measured  in  centimeters. 
Although  an  equation  like  (3.31)  is  cited  in  most  texts  as  the  reason 
oil  calms  troubled  waters,  the  enhanced  damping  is  really  not  very 
impressive,  especially  for  wavelengths  much  less  than  a  centimeter. 

To  examine  this  question  further,  we  have  solved  Equation  (3.27) 
numerically  for  a  range  of  different  values  of  B.  To  simplify  the 
discussion,  we  model  the  surfactant  by  an  ideal  gas,  and  imagine  that 
it  is  gradually  compressed  by  increasing  the  pressure  it  which 
confines  it.  Comparing  Equations  (2.1)  and  (2.2),  we  see  that  B  =  it 
in  this  approximation.  The  pressure  also  enters  the  surface  tension 
through  Equation  (2.1).  Experiments  which  measure  capillary  wave 

1  £_  1  Q 

damping  °  typically  use  a  transducer  at  a  fixed  frequency,  and 
measure  an  inverse  attention  length  k^  by  fitting  the  envelope  of  the 
wave  train  to  an  exponential  decay,  -  exp(-k(jx).  Numerical  solu¬ 
tions  of  Equation  (3.27)  at  a  fixed  real  wavevector  k  lead  to  a  com¬ 
plex  frequency  of  the  form 

to  =  to  (k )  -  iY ( k ) 
o 


(3.32) 


where  <joQ(k)  is  given  approximately  by  (3.28)  and  Y(k)  represents  the 

damping.  To  make  contact  with  the  experiments,  we  need  to 

analytically  continue  to  real  frequencies  and  complex  wavevectors. 

Let  k„  be  the  real  wavevector  such  that  u  (k  )  equals  the  real  fre- 
o  oo 

quency  w  in  question.  Replacing  k  by  kQ  +  ikd  in  Equation  (3-32), 
and  expanding  in  the  small  quantity  kd,  we  find  that 

kd  =  Y(k  )/[(du>  /dk)|  ]  .  (3.33) 

u  o 

The  group  velocity  rather  than  the  phase  velocity  enters  because  a 
spatially  narrow  transducer  at  fixed  frequency  generates  a  band  of 
wavelengths  in  the  vicinity  of  kQ. 

Figure  *4  shows  a  plot  of  kd  vs.  pressure  for  a  0.52  cm  capillary 
wave.  Also  shown  is  the  damping  for  an  incompressible  film.  The 
damping  increases  with  increasing  surfactant  concentration,  and 
remarkably,  goes  through  a  maximum  before  approaching  the  result  for 
an  incompressible  film  at  high  pressure.  This  nonmontonic  behavior 
has  been  observed  in  many  laboratory  experiments1^"1^,  including 

i  Q 

those  on  naturally  securing  sea  surface  films.  It  was  first 
pointed  out  theoretically  by  Dorrestein.^  Note  that  the  damping  can 
actually  exceed  the  result  for  an  incompressible  film.  As  noted  by 
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Davies  and  Vose  ,  this  must  be  due  to  an  appreciable  tendency  of  the 
horizontal  velocity  near  the  surface  to  reverse  relative  to  the  hori¬ 
zontal  velocity  in  the  bulk.  The  effect  of  a  finite  intrinsic 
surface  viscosity  is  to  reduce  the  maximum,  ultimately  leading  to  a 
monotonic  variation  in  the  limit  of  large  viscosities.^  Although  the 
enhanced  damping  at  intermediate  pressures  is  larger  than  that 
predicted  by  Equation  (3-30),  it  is  still  not  a  particularly  dramatic 


effect. 


H.O  EFFECT  OF  TURBULENT  FLUCTUATIONS  IN  THE  ATMOSPHERE 


The  linearized  hydrodynamic  theory  described  in  Section  3  is 
easily  generalized  to  allow  a  coupling  to  fluctuating  components  of 
the  Reynolds  stress  in  the  wind.  Upon  adding  the  zz  and  xz 
components  of  the  atmospheric  stress  tensor  p'  at  the  interface  to 
the  boundary  conditions  (3.13)  and  (3.23),  we  have 


('P  *  2n3zvz)|z=?  -a  g  =  p'z(x.t), 


(H.1) 


n(3  v  +  3  v  ]|  +  u)32v  =  p  (x,t), 

z  x  x  ZJ|Z=£  v OJ  x  x  xz 


(4.2) 


where  pzz’  (x,t)  and  pxz'  (x,t)  are  fluctuating  random  variables. 

Effects  due  to  the  mean  wind  velocity,  embodied  in  the  average  values 

of  pzz'  and  Pxz',  can  amplify  or  diminish  preexisting  capillary 

1  9 

waves,  depending  on  the  circumstances. 

Here,  we  study  generation  of  capillary  waves  due  to  f luctuat ions 
in  pzz'  and  pxz'  about  their  mean  values.  Assuming  a  wind  speed  of 
1  m/sec  at  a  height  of  1  meter  above  the  ocean,  we  estimate  that  tur¬ 
bulent  fluctuations  in  the  bulk  stress  tensor  should  be  present  down 
to  a  Kolmogorov  scale  of  about  .03  cm.,  which  includes  the  capillary 
wave  part  of  the  spectrum.  Because  little  is  known  about  such  fluc¬ 
tuations  close  to  an  air-water  interface,  we  shall  simply  assume  that 


the  fluctuating  components  of  the  wind  stress  are  unaffected  by  the 
presence  or  absence  of  surfactant,  and  compare  the  mean  square  wave 


amplitudes  generated  at  a  various  wave  vectors  and  frequencies. 


Consider  the  response  to  one  particular  Fourier  component  of  the 
fluctuating  wind  stress  tensor  at  the  interface, 


,  ,  ..  o  ikx-iuit  ,  ,  .  .  o  ikx-iu>t 

p'  (x,t)  =  p°  e  .  p'  (x,t)  =  p°  e 

Kzz  *zz  Kxz  ’  xz 


(4.3) 


where,  in  contrast  to  Section  3.0,  both  k  and  ui  are  real.  The 
analysis  of  Section  3.0  (again  setting  the  2d  viscosity  to  zero) 

now  leads  to  a  matrix  equation  for  the  amplitudes  A  and  C,  namely 


2  2  2  2  o 

~u)  +  to  "  2iYo)k  ,  iui  +  2\>u>k.l  ,  /A\  /-iup 

o  O  \  /  \  /  ZZ 


2vu)k2  +  iBk3/p  ,  u)2  +  2ivwk2  -  Bk2l/p'  '  c'  -iuip  °  .  (4.4) 

0  o  *xz 


The  amplitudes  which  follow  from  (4.4)  are 


A  =  -iw[(i(jj  2  +  2vukil)p  0  +  (2vwk2  -  Bk3i/p  )p  °]/A(k,w)  (4.5a) 

O  ZZ  O  X  z 


=  -  io)[  (u>2  «■  2ivwk2  -  Bk22,/p  )p  °+  (-u>2  +  u>  2)p  °]/A(k,w)  (4.5b) 

O  ZZ  o  Kxz  ’ 


where  A(k,w)  is  the  determinant  defined  in  Equation  (3.27).  We  can 
now  obtain  the  space  and  time  dependence  of  the  interface  height  via 


Equation  (3.9). 

Suppose  hat  Pzz°  and  Pxz°  are  chosen  from  an  ensemble  of  ran¬ 
domly  chosen  amplitudes,  and  assume  for  simplicity  that  pzz°  and  pxz 
ire  statistically  indepedent.  It  then  follows  from  Equations  (4.5) 
3.9)  that 

x.t:  =  3L(k,w)(p770)2  +  ST(k,u)(pxz°)2  (4.6) 

i-  aerates  an  ensemble  average  and 

.  *  :;v<!U  -  ia)2-  2vk2ui  -  iBk3i/pQ  |  2/  |  A  ( k ,  to)  |  2  (4.7a) 

3k”l/p0  *  iu)2_iu)02+  2vk2io|2/|A(k,u)  [2  (4  7b) 

3.  <,l'  are  response  functions  describing  how  the  inter- 

f  a- ■"  •  e  i  •*.  ,  v  •  •r.gi-.'idinal  '  i.e. ,  normal  to  the  surface)  and  trans- 

v.iveiy.  r'igure  5  shows  the  real  part  of  the 
r.  :  .H'  /  wave  iisper3ion  relation  (3.26).  It  is 


instructive  to  imagine  fixing  the  wavevector  and  varying  frequency 
along  the  vertical  line  shown  in  the  figure.  The  response  functions 
and  ST  should  become  large  when  the  frequency  passes  through 
u)  (k).  The  width  of  this  resonance  is  proportional  to  the  damping 
discussed  in  Section  3-0.  The  height  of  the  resonance  is  a  measure 
of  how  easy  it  is  to  generate  waves  at  the  resonance  frequency. 

Figures  6a  and  b  contrast  the  response  with  and  without  a 
surface  film  in  the  gravity  wave  part  of  the  spectrum.  The  wave- 
'ength  is  one  meter,  we  have  assumed  that  the  surface  film  has  a 
compressibility  of  10  dyne/cm.  On  a  clean  surface,  there  is  an 
extremely  narrow  resonance  at  a =  7.85216  Hz.  The  surfactant  shifts 
the  resonant  frequency  and  leads  to  a  small  increase  in  the  damping, 
as  measured  by  the  half-width  at  half-maximum.  The  heights  of  the 
resonances  are  depressed  only  slightly  by  the  surface  film.  Figure  6 
shows  that  there  will  be  persistent  long  wavelength  "oily  swells"  in 
sea  slicks,  notwithstanding  the  presence  of  a  surface  film.  This 
result  is  to  be  expected,  since  the  right-hand  side  of 
Equation  (3.23)  becomes  negligable  compared  to  the  left-hand  side  at 
sufficiently  long  wavelengths. 

The  results  for  capillary  waves  with  wavelength  0.52  cm  are 
shown  in  Figures  7a  and  7b.  Here,  the  effect  of  a  surface  film  is 
much  more  remarkable.  Particularly  striking  is  the  reduction  in  the 
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Figure  6a.  Longitudinal  response  functions  in  the  gravity  wave  part  of  the  spectrum 
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Figure  6b.  Transverse  response  functions  in  the  gravity  wave  part  of  the  spectrum 
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Figure  7a.  Longitudinal  response  functions  in  the  capillary  wave  part  of  the  spectrum. 


Figure  7b.  Transverse  response  functions  in  the  capillary  wave  part  of  the  spectrum 


heights  of  the  response  functions.  The  longitudinal  response  is  down 
by  a  factor  of  20  and  the  transverse  response  is  barely  visible  (down 
by  a  factor  of  60).  A  significant  increase  in  the  damping  is  also 
apparent,  as  well  as  a  small  shift  in  the  resonant  frequency 
(c =  369.06  Hz  for  a  clean  surface).  Figure  7  shows  that  oil  calms 
troubled  waters  not  only  because  of  this  increased  damping,  but  also 
because  it  is  very  difficult  to  generate  large  amplitude  capillary 
waves. 


The  reduced  amplitude  for  capillary  waves  is  most  pronounced  in 
the  transverse  response  function.  In  contrast  to  the  damping,  which 
approaches  a  finite  limit  for  an  incompressible  film,  the  amplitudes 
generated  by  transverse  fluctuations  actually  vanish  as  B  tends  to 
infinity.  To  see  this  in  more  detail,  we  start  with  the  time- 
dependent  boundary  conditions  (3.13)  and  (3.23),  including  the  effect 
of  a  transverse  wind  stress  pxz(x,t), 


2n  3Z  vx  ♦  p  gc  ♦  P0  |f  -  0  M  =  °’ 

3x 


(4.8a) 


n(3  v  ♦  3  v  )  -  —  3  v  =  p  (x,t), 
x  z  z  x  0J  xx  Kxz 
0 


(4.8b) 


where  we  have  set  u2  =  0  and  used  Equation  (3.8)  for  the  pressure. 
All  quantities  in  Equations  (4.8)  are  to  be  evaluated  with  z  =  5. 
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The  effect  of  viscosity  in  these  boundary  conditions  can  be  neglected 
in  the  large  B  limit.  Upon  setting  n  =  0,  using  (3.1)  and  differen¬ 
tiating  each  equation  with  respect  to  time,  we  find  Equations 
for  4>  and  <Jj,  namely 


P0gOz4>  "  3X<I>)  +  Pq3£  4>  -  o  3x2  Oz*  -  3^)  =  0  (4.9a) 


B  3  2  (3  <t>  +  3  41)  =  3.  p 
x  x  zr  t  *xz 


(4.9b) 


We  now  assume  that  41  and  \Jj  can  be  written 


♦  -  A ( t )  eikx  +  kZ  -  iuo(k) 


(4.10a) 


.  r.  .  ikx  +  Hz  -  iw  (k) 
4/  =  C(t)  e  o 


(4.10b) 


where  k  is  real,  w  (k)  is  given  by  Equation  (3.26),  and  A(t)  and  C(t) 
are  slowly  varying  functions  of  time.  The  wave  represented  by 
Equation  (4.10)  will  be  acted  upon  most  strongly  by  that  Fourier 
component  of  pv„  whose  wavevector  and  frequency  match  k  and  co  (k). 

Kl.  O 

1  Q 

As  discussed,  for  example,  by  Longuet-Higgins  ,  the  wave  will  be 
most  strongly  amplified  when  the  phase  of  pxz  makes  this  variable 
stress  a  maximum  at  the  wave  crests  and  a  minimum  at  the  troughs. 
Assume  for  simplicity  that  only  this  mode  is  present,  and  that  it 
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varies  as  in  Equation  (4.3).  When  Equations  (4.10a)  and  (4.10b)  are 
substituted  in  (4.9),  the  result  is  two  ordinary  differential  equa¬ 
tions  for  A(t)  and  C(t),  namely 


A(t)  +  iA(t)/2u>  +  ^  oi  C(t)  =  0 

o  2  o 

ik^A(t)  +  k2HC(t)  =  -iio  p  °. 

o'xz 


(4.11a) 

(4.11b) 


We  now  confine  our  attention  to  the  capillary  wave  part  of  the  spec¬ 
trum.  Neglecting  the  second  order  time  derivative  in  (4.11a),  and 
using  (4.11b)  to  eliminate  C(t),  we  find  a  simple  differential  equa¬ 
tion  for  A(t)  (ignoring  a  contribution  to  the  imaginary  part  of  A), 
namely. 


dA 

dt 


J _ 

2/*2 


2  % 


(^L) 
Po 


,  7/4  . 
k  A  + 


1  (ok^3/4  /J 

2/~2  Po  Bk2 


(4.12) 


If  Pxz°  =  0.  the  wave  dies  out,  and  A ( t )  decays  with  the  damp¬ 
ing  rate  of  an  incompressible  surfactant,  given  by  the  imaginary  part 
of  Equation  (3*30).  For  a  nonzero  tangential  stress,  A(t)  approaches 


.  surf 
A 


p  °/Bk3 
*xz 


(4.13) 


for  long  times,  a  limit  which  tend3  to  zero  for  large  B. 


Equation  (4.12)  should  be  contrasted  with  the  result  for  a  clean 
1 Q 

surface,  which  is 


^  -  2vk2A  +  p  °/2p  . 

dt  xz  o 


(4.14) 


This  equation  has  the  long  time  limit 


.clean  o,,.  ,2 

A  =  p  /4vk  p  . 

®  *xz  o 


(4.15) 


The  ratio  Aclean/Asurf  can  t>e  made  arbitrarily  large  as  B  becomes 

00  CO 

large,  i.e.,  for  sufficiently  incompressible  films. 
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